Stability of fractional vortex states in a two-band mesoscopic superconductor 
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We investigate the stabihty of noncomposite fractional vortex states in a mesoscopic two-band 
superconductor within the two- component Ginzburg-Landau model. Our analysis explicitly takes 
into account the relationship between the model parameters and microscopic material parameters, 
such as partial density of states, fermi velocities and elements of the electron-phonon coupling 
matrix. We have found that states with different phase winding number in each band (Li ^ L2) 
and fractional flux can exist in many different configurations, including rather unconventional ones 
where the dominating band carries larger winding number and states where |Li— L2I > 1. We present 
a detailed analysis of the stability of the observed vortex structures with respect to changing the 
microscopic parameters, showing that, in the weak coupling case, fractional vortex states can be 
assessed in essentially the whole range of temperatures and applied magnetic fields in which both 
bands are active. Finally, we propose an efficient way of increasing the range of parameters for which 
these fractional vortex states can be stabilized. In particular, our proposal allows for observation of 
fractional vortex structures in materials with stronger coupling, where those states are forbidden at 
a homogeneous field. This is accomplished with the help of the stray fields of a suitably prepared 
magnetic dot placed nearby the superconducting disk. 



I. INTRODUCTION 

Many properties of the electronic condensate in a 
macroscopic quantum system can be unveiled by rotat- 
ing it and studying the structure of the induced quantized 
vortices i In particular, the amount of magnetic flux car- 
ried by Cooper-pair vortices in bulk conventional super- 
conductors determines whether the material is of type 
I or type 11.-^ In the first case, vortices tend to merge 
into large flux domains carrying many flux quanta. Con- 
versely, for the type-II material, vortices are singly quan- 
tized and repel each other, tending to arrange themselves 
into a triangular lattice. 

Some superconductors, however, are provided with two 
or more electronic condensates arising from Cooper pair- 
ing in different bands of the material. These multiband 
systems have attracted much interest in the last decade 
because: (i) most of the recently discovered supercon- 
ducting materials, such as MgB2 and compounds of the 
iron-pnictide family, are recognized as multiband super- 
conductors; and (ii) they exhibit a variety of new and 
interesting phenomena with no counterpart in conven- 
tional single-component superconductors. Perhaps the 
most intriguing of these phenomena are related to the 
exotic vortex structures that can emerge in a multicom- 
ponent superconductor. 

In the bulk, a multicomponent vortex can exist in equi- 
librium only in a so-called composite state where vortices 
in the different condensates share the same core.^ How- 
ever, because of the different length scales at which 
the Cooper-pair density varies in each component, the 
interaction between such composite vortices can be non- 
monotonic: short-range repulsive with an attractive tail^^ 
This ultimately leads to vortex clustering at low magnetic 



fields and, thereby, to the formation of a semi-Meissner 
state, that is, a mixture of flux- free and vortex-cluster re- 
gions coexisting in equilibrium. This phase was predicted 
by Babaev and Speight and experimentally observed by 
Moshchalkov et al& in an ultraclean MgB2 single crys- 
tal. In this case, the first critical field corresponds to the 
thermodynamic stabilization of a vortex cluster in the 
sample, rather than a single vortex. In this regime, be- 
sides clusters, other exotic configurations are likely, such 
as vortex rings, giant-multi vortex groupings, and other 
unusual patterns. 

Actually, the competing interactions in a multiband 
material can lead to even more unconventional vortex 
topologies. For instance, "unmatched" composite vor- 
tices, that is, with both condensates having different 
phase winding numbers, can exist out of equilibrium . ^^^^^ 
Such vortices enclose arbitrary magnetic flux, in contrast 
to quantized Abrikosov vortices in conventional super- 
conductors. Another example are noncomposite vortices, 
in which the phase singularities in each condensate are 
displaced from one another. These vortices also carry 
arbitrary flux, however, since their energy increases pro- 
hibitively with sample size, they can only be observed 
in small samples. These exotic vortex structures could 
be used as a hallmark of multiband superconductivity. 
Therefore, an important question is how to stabilize and 
detect such fractional vortices. 

Recently, Chibotaru and co-workers ^^^^^ have proposed 
that fractional vortices can be realized and even ther- 
modynamically stabilized in a two-band mesocopic disk. 
Subsequently, Geurts et al}^ extended the analysis by in- 
cluding magnetic coupling, besides the interband Joseph- 
son coupling used in previous works. All these stud- 
ies were performed within the two-component Ginzburg- 
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Landau (TCGL) formalism, taking phenomenological 
constants as parameters of the simulations. However 
as recently observed in Ref. 15, the phenomenological 
constants of the TCGL model cannot be chosen freely 
because they are coupled by microscopic material pa- 
rameters, such as electron-phonon coupling constants, 
Fermi velocities and partial density of states. A more 
systematic study of the stability of composite and non- 
composite vortex states in mesoscopic two-component 
superconductors, within the correct microscopic frame- 
work, is still pending. 

In this work, we investigate the stability regions of frac- 
tional flux vortices in a two-band mesoscopic supercon- 
ductor in a parameter space defined by microscopic ma- 
terial parameters. In addition, we demonstrate how the 
stability of these states can be considerably enhanced by 
the presence of a closeby magnetic dot with a suitable 
magnetization. 



II. THEORETICAL FORMALISM AND 
SIMULATION DETAILS 

We consider a mesoscopic two-band superconducting 
disk of thickness d much smaller than the penetration 
depth A and the two characteristic lengths of density vari- 
ations in both condensates, in such a way that the system 
is effectively two-dimensional. Our calculations rely on 
the minimal two-component Ginzburg-Landau (TCGL) 
model, where the only coupling appearing explicitly in 
the free-energy functional is Josephson-like.^^^^ Within 
this framework, the total free energy can be written as 
the sum of the free energies of the otherwise isolated 
bands and the Josephson coupling free energy, that is. 



2hujDe ^ ^ where a;/) is the Debye frequency and 



F 



+a,\<i>j\^ + -l3j\^j\ 



2e 

-ihV - —A 1 



r(^i^* + ^t^2) > (1) 



j = 1 , 2 is the band label and rrij , aj , j3j , and F are pa- 
rameters derivable from the microscopic theory. The last 
term accounts for magnetic coupling between the bands, 
where Hq is the applied magnetic field, and h = rotA 
is the local (total) field. However, because of the small 
thickness and lateral size of the system under investiga- 
tion, screening effects are negligible and thereby so is the 
magnetic coupHng between the bands. 

Within the two-band Eilenberger formalism,— i^J^ the 
relevant microscopic parameters are the Fermi velocities 
Vj^ the partial densities of states nj7V(0), and the ele- 
ments of the electron-phonon coupling matrix 
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are the roots of the linear system of self-consistecy equa- 
tions for the gaps. The correct value of Tc is given by the 
smallest root. It is also convenient to define the follow- 
ing auxiliary parameters: r = — In (T/Tc) 1 — T/Tc, 
W = 2\/27rTc/ Y^7C(3), r] = det A/nin2, and the positive 
constants 

Si = A22 — nir]S and S2 = An — n2r]S. 

In terms of the microscopic parameters above, 
the TCGL coefficients can be expressed as: aj = 



-N{0)nj{r-Sj/njrj), 



N{0)nj/W^, ruj 



3W^/{N{0)njvj), and F = N{0)Xi2/r]. As usual, one can 

define the healing lengths = hvj / V^W, which are re- 
lated to the characteristic lengths of density variations of 
order parameters 1 and 2 (see for instance Ref. 4). How- 
ever, the real healing lengths of the condensates have to 
be calculated, and are strongly influenced by the crit- 
ical temperatures of the otherwise uncoupled bands^^- 
(that is, the temperature at which the corresponding aj 
changes sign). 
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Notice that both Tcj are always smaller than Tc. There- 
fore, for Tci,Tc2 < T < Tc, both aj are positive and 
superconductivity survives in the system only due to cou- 
pling between the bands, whereas for T < Tci,Tc2, both 
bands are active {aj < 0). For temperatures such that 
TcP <T < TcA, one of the bands (band A) is active while 
the other (band P) is passive, that is, it remains super- 
conducting only because of Cooper pairs coming from 
band A. 

Finally, minimization of the TCGL free energy [1] leads 
to the following dimensionless GL equations: 

{-iV - Af^r - (n - |V^ip)V^i - = 0, (3a) 

niT] 

4(-^V - AflP2 - {t2 - \^2?)^2 - —^1 = 0, (3b) 

v{ n2r] 

-ln{T/Tcj) = r - Sj/ifijf]). Here, we 
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The critical temperature is given by the relation 1.76Tc = 



adopted the following temperature independent units: Tc 
for temperatures, W for both order parameters, ^1 for 
distances, and Aq = he/ 2e£^i for the vector potential. 
Hereafter, we choose the disk radius R = lO^i <C Aeff, 
where Aeff = A^/d is the effective penetration depth, 
so that screening and demagnetizing effects can be ne- 
glected. Accordingly, the vector potential in (|3a|) can be 
well approximated by that corresponding to the external 
magnetic field. 

The minimal TCGL model with the r 1 — T/Tc 
approximation is in principle strictly valid only in the 
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immediate vicinity of Tc, in which case the model is re- 
duced to the conventional GL theory featuring a single 
order parameter with a single length scale. ^- However, 
for lower temperatures the correct microscopic scenario 
points to two condensates with two a priori different 
healing lengths. Inspired by the success of the GL for- 
malism in producing correct results well out of its strict 
validity range in several other instances, one may ask 
whether the applicability range of the TCGL model can 
be pushed to lower temperatures. Indeed, in Ref. JJi it 
was shown that an extension of the GL theory to include 
complete r^/^ terms to the expansion of the order param- 
eters pushes the validity of the theory to temperatures as 
low as 0.62Tc, depending on the values of the microscopic 
parameters. In addition, this model correctly predicts 
two different density length scales, which become equal 
only for T ^ Tc, thus in reconciliation with previous 
results from the microscopic theory. Unfortunately, the 
extended GL model for the non-zero field case is to date 
available only for the single-band superconductors.^^ 

Recently, evidences that the standard TCGL model 
provides a correct description of two-band superconduc- 
tors has been given in Ref. H The authors set out a 
systematic comparison with the microscopic two-band 
Eilenberger theory suggesting that the minimal TCGL 
model can offer a quite accurate quantitative description 
of features specific of two-band superconductivity, includ- 
ing two- length-scales vortex solutions, if one retains the 
full logarithmic temperature dependence of r for temper- 
atures smaller but not too close to Tc. Good agreement 
was achieved for temperatures as low as T = 0.85. There- 
fore, to ensure the validity of the TCGL model used in 
the present work, we fix the temperature of the sample at 
T = 0.90 and choose microscopic parameters such that 
both Tci and Tc2 are also close to Tc (typically, higher 
than 0.9). Such critical temperatures are also chosen in 
order to minimize the effects of hidden criticality in the 
case of weak coupling between the bands (see Ref.^^ for 
details). 

The calculations are performed as follows. Eqs. [3alare 
numerically integrated by a relaxation method where a 
diffusion- like time-dependence of the fields is assumed. 
This procedure is similar to solving the time-dependent 
Ginzburg-Landau equations. However, once we are in- 
terested only in the stationary states, we chose the same 
diffusion constant, D = 1, for both order parameters. 
The equations are then discretized in space following the 
link- variable method and the resulting gauge- invariant, 
finite-difference equations are time-integrated via a semi- 
implicit scheme. As compared to the conventional Euler 
method, the semi-implicit integration provides better sta- 
bility and higher accuracy, which are particularly crucial 
for the convergence of fractional vortex solutions. The 
results are presented as functions of the external mag- 
netic flux threading the disk area, ^ = H/ttR^, in units 
of the magnetic flux quantum, = hc/2e. The config- 
urations corresponding to bands 1 and 2 having a phase 
winding number Li and L2, respectively, are labeled as 



III. INFLUENCE OF MICROSCOPIC 
PARAMETERS 

In most previous studies on fractional vortex states, 
the main focus of analysis was on specific material pa- 
rameters. However, to date there is no consensus about 
the correct coupHng constants of, for instance, MgB2, 
by far the most studied two-band material (see for in- 
stance Ref. 22 and references therein). Therefore, here 
we adopt a different strategy and study the possible con- 
figurations of Cooper pair densities in different regions 
of the parameter space defined by the microscopic quan- 
tities A12, ni (= 1 — 712) and V2. Notice that in our 
calculations vi is a constant bound to the unit length. 
For a treatable analysis, we fix only the intraband cou- 
pling constants. An = 2.415 and A22 = 1.211. The only 
bias in such a choice is that these values lead to Td and 
Tc2 close to Tc, which is important for the applicability 
of the Ginzburg-Landau framework. However, we believe 
that other values would impart no qualitative change to 
the main conclusions of our investigation. 



A. Effect of the partial densities of states 

Here, we fix the interband coupling at A12 = 0.001 and 
the Fermi velocity ratio at Vi/v2 = 1.225 and change 
only the partial densities of states, ni and n2. The main 
effect of changing these quantities is the relative vari- 
ation of the critical temperatures of band activity, Td 
and Tc2 • Although our choice of intraband coupling con- 
stants (with All — 2A22) benefits band 1 as, in principle, 
the strongest band, one can always have Tc2 > Td if ni 
and n2 are properly chosen. In fact, it is the product 
UjXjj which determines which of the bands will domi- 
nate. It results from Eq. [2]that both bands has the same 
Tcj when niAn = n2\22- For our choice of the coupling 
matrix, this reads ni = 1 — n2 = 0.334. Band 1 is the 
dominating band above this line only, whereas band 2 
dominates below it. By moving away from the Td = Tc2 
Hue, in any direction, one of the bands will become pas- 
sive. This happens with band 1 for ni = 1 — n2 < 0.317 
and with band 2 for ni = 1 — n2 > 0.354. 

Fig. [1] presents the stability regions of (^1,^2) states 
in the ni-^ plane. The curves plotted in these diagrams 
correspond to the limit of stability of each particular 
(Ti,T2) state, as explained in the caption. For a given 
ni, each stationary state was obtained after initializing 
the system with a particular (^1,^2) configuration at a 
given value of $ and checking whether it converges to 
a stationary solution. Once a field value was found in 
which that configuration is stable, we swept the field up 
and down to check the whole range of stability. As ex- 
pected, the appearance of fractional states, i.e., states 
with Li 7^ T2, is only possible in the range of ni and 
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FIG. 1. (Color online.) Stability regions of integer and frac- 
tional (Li, L2) configurations in the plane defined by the par- 
tial density of states m and reduced external magnetic field 
$/$o. Integer (L, L) states are bounded by two full lines, the 
left (right) one corresponding to the smallest (highest) field 
value below (above) which the state becomes unstable. The 
stability regions of fractional states correspond to the closed 
regions delimited by full lines for Li < L2 and dashed lines 
for Li > L2. For comparison, we show the stability limits 
of integer states as light gray lines. The heavy gray lines in 
the upper panel delimits the region (shaded area in the lower 
panel) where both bands are active. 



712 values for which both bands are active. Outside this 
region, the active band forces the passive band to follow 
its configuration via Josephson coupling. Surprisingly, 
both combinations Li < L2 and Li > L2 appear as sta- 
ble (1/1,^2) states in most part of the two- active- band 
region, no matter which of the two bands has a higher 
Tci. This is in contrast with previous studies, where only 
states with Li < L2 (band 1 being the strongest band) 
were reported. 

A detailed analysis of how these states change with 
field and how they are energetically compared to one an- 
other is presented in Fig. [2l The results shown in the 
upper panel of this figure were calculated for ni = 0.343, 
for which Td = 1.0 and Tc2 = 0.95. Therefore, the "in- 
verted" states are those with Li > L2. The variety of 
possible configurations in a two-band material is imme- 
diately revealed. For instance, at a field ^ = 3.2<l>o we 
found as much as seven different stable configurations. 
It is also clear in this figure that all fractional states 
are metastable, i.e. have higher energy than the ground 
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FIG. 2. (Color online.) Free energy as a function of external 
magnetic field (flux) . Full lines correspond to states stabilized 
by sweeping the field up and down starting from the Meissner 
state at if = as initial condition (see text). Dashed lines 
correspond to the states obtained only by imposing a suitable 
initial condition, as explained in the text. Top: all states 
found for the case n\ — 0.343. Bottom: selected fractional 
states [(0,1), (1,0), (0,2) and (2,0)] for three different values of 
the partial density of states: n\ — 0.343 (left), 0.334 (middle), 
and 0.325 (right). 



state. Note that fractional states can be also found in the 
ground state, but for very weak interband coupling, as 
discussed later in Sec. lIII Cl In the present case, the states 
with L\ < L2 have, in general, lower energy, sometimes 
very close to the ground states. An exception is the (0,2) 
state, with a remarkably high energy. The evolution of 
relative energies and stability range of a selection of frac- 
tional states as one decreases ni is illustrated in the lower 
panel of Fig. El As expected, for m = 0.325 (T^i = 0.95 
and Tc2 = 1-0) the roles of states with Li > L2 and those 
with Li < L2 are interchanged, with the former having 
in general lower energy. 

Another important feature distinguishing "regular" 
and "inverted" (^1,^2) solutions concerns the way these 
states are found. The curves in Fig. [2] represented by 
full lines were obtained after initializing the system in 
either the Meissner state, at if = 0, or normal state, 
at H just above the sample's upper critical field, and 
then isothermally sweeping the field up and down. Such 
a procedure accounts for a realistic magnetic history 
and allows one to access both metastable and ground 
states. Experimentally, it is similar to performing minor 
isothermal magnetization-loops measuremets on meso- 
scopic samples^^^^^. In this case, finding a new vor- 
tex state involves entry or exit of one or more vor- 
tices through the surface barrier once a saddle point 
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is reached Therefore, we anticipate that the frac- 
tional states represented by fuh hnes in Fig. [21 though 
metastable, should be easily accessed experimentally. On 
the other hand, those states represented by dashed lines 
could not be obtained this way. Instead, we used a par- 
ticular vortex configuration as the initial state, in an- 
ticipation of a similar stable outcome of the calculation. 
We are not aware of an experimental counterpart of this 
process. One could then conclude that, although in- 
verted fractional states do appear as stable solutions of 
the Ginzburg-Landau equations, they are unlikely to be 
observe experimentally. 



B. Effect of the Fermi velocities 

In this subsection, we retain the coupling constants 
defined in the previously, fix ni = 1 — n2 = 0.343, and 
analyze the effect of Fermi velocities on the vortex states 
in a two-band mesoscopic disk. As shown in the theoret- 
ical formalism, the Fermi velocity of a band is directly 
proportional to the healing length of the corresponding 
condensate at zero temperature. In other words, Vj is 
closely connected to the maximum critical field band j 
can sustain in the zero coupling limit and at T = 0. 
Bearing in mind that is fixed as our unit length (and 
so vi remains fixed as well), by changing the ratio Vi/v2 
as a parameter we are actually changing ^2- Therefore 
the critical field above which band 2 becomes passive is 
expected to increase with {vi/v2)'^. Fig. [3] presents the 
evolution of the stability regions of fractional states as a 
function of the squared Fermi velocity ratio {vi/v2)'^. In 
addition, we plot the phase boundary separating the ac- 
tive and passive states of band 2. Above this line, both 
bands are active. The critical field for band 1 activity 
(not shown) corresponds to the constant value 9.5<l>o. 

The role played by the Fermi velocities and the healing 
lengths in determining the shape of the stability regions 
in Fig. [3] can be identified as follows. As one approaches 
the phase boundary for band 2 activity, ^1 increases lin- 
early and so does the size of a vortex in that band. This 
is illustrated in Fig. [3] (b) and (c), where we plot the 
Cooper pair distributions of band 1 in the (1,0) state 
and band 2 in the (0,1) state at a fixed applied field 
and different values of {vi/v2)'^. As an estimate, we de- 
fine the vortex size as the contour where the absolute 
value of the order parameter recovers by 50% its maxi- 
mum value in the zero coupling limit, which is given by 

1 /2 

tpjo = Wt- . Within this criterion, one can clearly ob- 
serve the expansion of the band-2 vortex in the (0,1) state 
as one approaches the phase boundary. Such an increase 
in the characteristic length of density variations in band 
2 induces a decrease in the energy barriers in this band. 
Hence, when approaching this phase boundary, all pro- 
cesses involving entry or exit of a vortex in band 2 should 
take place in advance (at a lower field for increasing ^ 
and higher field for decreasing On the other hand, 
since is fixed, one would expect no change in the size 
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FIG. 3. (Golor online) (a) Stability regions of fractional 
(Li,L2) states in the plane defined by the ratio between 
the Fermi velocities {vi/v2) and the applied magnetic flux 
The gray line corresponds to superconducting-normal phase 
boundary of band 2 in the absence of interband coupling, (b) 
and (c): Plots of the absolute value of the order parameter 
of band 1 in the (1,0) state (b) and band 2 in the (0,1) state 
(c) for the point I, II, and III indicated in (a). The order pa- 
rameters are normalized by the maximum zero-coupling value 
of each band. The circles indicate the size of the vortex as 
deflned in the text. 



of a vortex in band 1 [which is indeed approximately the 
case as observed in Fig. [3] (c)] and thereby no important 
change in the energy barriers of that band. In this case, 
processes involving entry or exit of vortices in band 1 are 
expected to take place at field values that do not change 
considerably with the parameter {vi/v2)'^' 

The scenario explained above has an immediate conse- 
quence on how the supercooling, ^sci and superheating, 
^sh^ fields of each fractional state changes with {vi/v2)'^. 
In general, by sweeping the field up and down, both 
(L, L+1) and (i^+1, L) states will decay to (i^+1, L^l) at 
the corresponding superheating field in the upward sweep 
or to (L, L) at the corresponding supercooling field in the 
downward sweep. When the original state is {L^L + 1), 
the (L + 1,L + 1) state is achieved with the entry of a 
vortex in band 1 and is essentially constant, while 
(L, L) is obtained via the escape of a vortex from band 2, 
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FIG. 4. (Color online.) The stability regions of fractional 
(Li,L2 / Li) states in the plane defined by the inter band 
coupling coefficient A 12 and the magnetic flux $ through the 
sample area. In shaded regions fractional states represent the 
thermodynamic equilibrium of the studied system. 



thus pushing ^sc to higher values as approaching band-2 
phase boundary. Conversely, when (L + 1, L) is the orig- 
inal state, ^sh is pushed to lower values, because of the 
advance entry of a vortex in band 2, while ^sc keeps es- 
sentially constant, once only the energy barrier of band 
1 is probed. 



C. Effect of the interband coupling 

The effect of Josephson coupling on the stability of 
fractional states has been studied in Ref. and, in more 
detail, in Ref. J4. In these works the Josephson cou- 
pling strength T was treated phenomenologically, with- 
out taking into account that the microscopic, interband 
coupling constant A12 can also exert an important influ- 
ence on Tcj and Tc, specially for higher coupling values. 
In this subection, we investigate the stability regions of 
fractional states in the Ai2-^ plane in consonance with 
the microscopic background and extend the analysis to 
the cases of (L + configurations. Our main results 
are featured in Fig [H The main trends shown in this 
figure are similar to those presented in Ref. shrink- 
ing of the stability regions as A12 increases and frac- 
tional states of lower vorticity are capable of surviving to 
stronger coupling. Inverted fractional states also present 
similar behavior. However, their regions of stability are 
displaced to lower field values. This can be understood 
in light of our explanation of the supercooling and su- 
perheating processes of {L + l^L) and {L^L + 1) states 
(c.f. Sec. line]) . Here, band 2 is always the weaker band. 
Therefore, ^sh{L + 1,1/) should be advanced to a lower 
value because it involves vortex entry in band 2, while 
^sc{L + 1,L) should be delayed to a lower value since it 
involves escape of a vortex from the stronger band 1. 



Another particularity of (L + 1,L) states is that, be- 
cause of their higher energy, they are less resilient than 
(L, L + 1) states. In fact, for the set of parameters consid- 
ered here, (L + 1, L) states can never be found in thermo- 
dynamic equilibrium. In contrast, regions of thermody- 
namic equilibrium for (L, L + 1) can be observed at very 
low values of A12, as can be observed in Fig. IH where 
these regions are represented by shaded areas. 



IV. INFLUENCE OF AN OFF-PLANE 
MAGNETIC DOT 

In this section, we analyze the effects of an inhomo- 
geneous field on the vortex structures observed in the 
previous section. We start with particular microscopic 
parameters for the superconducting disk (An = 2.415, 
A22 = 1.211, A12 = 0.001, V1/V2 = 1.225 and m = 
1 — n2 = 0.343) such that fractional states are acces- 
sible. The inhomogeneous field is provided by a cylin- 
drical magnetic dot, placed coaxially with the supercon- 
ducting disk and having a permanent, homogeneous off- 
plane magnetization given by M = z^m/t^R'^- (Here, 
is not to be confounded with the flux generated by 
the magnetic dot through the superconductor. Rather, it 
is only a convenient way to express magnetization in the 
same units as magnetic field.) To simplify our analysis 
we chose geometric parameters such that the stray fields 
generated by the magnetic dot at the superconducting 
disk have a flat distribution. This is accomplished by a 
magnetic cylinder of radius Rmd = O.bR = 5, height 
h = 4.8, and placed a distance / = 0.7 above the super- 
conducting disk (with assumed oxide layer in between, to 
avoid proximity effect). 

Fig. [5] illustrates the evolution of the stability regions 
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FIG. 5. (Color online.) $m-$ diagram showing the stability 
regions of fractional (Li, L2 ^ Li) states in the plane defined 
by the magnetization of the magnetic dot and the external 
magnetic field, both multiplied by the superconductor area. 




(a) (b) (c) (d) 

FIG. 6. (Color online.) Absolute value of the order parameters and external field profiles (provided by the magnetic dot and 
homogeneous field source) for the (2,3) state at the median of its stability region for (a) $m = 2.25 and $ = 2.0, (b) $m = 
and $ = 5.15, (c) $m = —0.75 and $ = 6.05, and (d) $m = —1.125 and $ = 6.33. The order parameter of each band is 
normalized by its respective zero-coupling maximum value and plotted in log scale for better visualization. Dashed circles 
represent the edge of the magnetic dot. 



of (L, L + 1) and (L + 1, L) states (for < L < 4) when 
the magnetization of the disk is changed. It is promptly 
observed that the stability of all fractional states can be 
considerably enhanced by the presence of the magnetic 
dot. This happens for moderate positive values of ^m- 
Too large positive values tend to drive the entire sys- 
tem to the normal state, which of course also destroys 
the fractional states. On the other hand, negative is 
deleterious to fractional states because in that case the 
central region of the superconducting disk encloses less 
flux, which disables it to host vortices in that area - while 
disk periphery is in any case not favorable for vortices to 
reside. Therefore, as the external field in this region is 
fully compensated, the Meissner state becomes predomi- 
nant. Upon further decreasing to even more negative 
values, the system becomes susceptible to the appearance 
of integer and fractional antivortex states. The first frac- 
tional states appearing are (0,-1) and (—1,0), followed 
by (-1,-2) and (-2,-1), and so on. The stability re- 
gions of such fractional antivortex states correspond to 
a rotation of the stability regions of their counterparts 
with positive vorticity, shown in Fig. [5l by 180° with 
respect to the origin. For the parameters used in our 
simulations, we did not detect any states where vortices 
and antivortices coexist in neither the same band nor in 
different bands, though they are expected to appear for 
larger samples as already demonstrated in single band 
superconductors 

The role of magnetic moment orientation can be bet- 



ter understood by analyzing the field profiles for different 
$M and $ and the corresponding vortex configurations. 
For such an analysis, we focus on a particular fractional 
state, (2,3), and study the evolution of the Cooper-pair 
density when $m and <l> are simultaneously changed in a 
way as to keep the (2,3) state at the median of its stabil- 
ity field range for each magnetization value. The results 
are presented in Fig. [6] where field profiles and Cooper- 
pair distributions of both bands are shown for the cases 
^m/^o = 2.25, 0, -0.75, and -1.125. In all cases the lo- 
cal flux in the central region of the superconducting disk, 
where the vortices sit, has approximately the same value, 
~ 5<l>o, which is roughly the flux necessary to stabihze 
state (2,3) for <I>m = 0. Therefore, the energy barrier for 
vortex escape is essentially the same in all cases in such 
a way that the supercooling field is expected to be an 
approximately linearly decreasing function of On 
the other hand, for a new vortex to come in, it has to 
probe the energy barrier near the superconductor edge, 
where the flux profile is very sensible to the magnetiza- 
tion of the dot. For positive ^m, screening currents are 
weaker near the sample edges making it more difficult 
for a new vortex to enter. Therefore the energy barrier 
for vortex entrance increases (with respect to the median 
of the stability range) with Hence, superheating is 
prolonged to higher field values, thereby providing the 
broadening of the stability range of a given vortex con- 
figuration. Note that this picture holds for any vortex 
arrangements, including integer states. 
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FIG. 7. (Color online.) Effect of the magnetic dot magne- 
tization on the stability of state (0,1) in the Ai2-$ plane for 
different values of the reduced magnetization $m. 



A striking feature that can be observed in Fig. [6] is the 
expansion of vortex configurations when <I>m is decreased 
down to negative values. Expansion of very compact vor- 
tex configurations with the help of ferromagnetic dots 
has been anticipated to occur for single-band mesoscopic 
squares in Ref. |3l| In the present case, , we demonstrate 
that this effect can also be applied to multiband mate- 
rials, thus providing an invaluable tool for experimental 
visualization of fractional vortex states. This expansion 
occurs because the decrease of the magnetic flux in the 
central region of the disk reduces the confining power of 
the screening currents and allows vortices to repel each 
other further from the center of the sample. Further re- 
duction of will lead to the sequential expulsion of 
vortices until the Meissner state is reached. 

In addition to discussing the influence of the magnetic 
dot on the stability regions for given interband coupling, 
it is also important to analyze what happens when A12 
is varied. Fig. [7] demonstrates that not only the range 
of stability of fractional states can be considerably en- 
hanced but the maximum value of A12 in which a given 
fractional state is stable increases as well. That is to 
say, a properly magnetized magnetic dot placed close to a 
mesoscopic two-band superconductor can optimally stabi- 
lize fractional vortex states, allowing them to be observed 
even in stronger intraband-coupling materials where they 
would be otherwise forbidden in a homogeneous field. For 
the (0, 1) case illustrated in Fig. Owe found that the op- 
timum magnetization is <I>m ^ 1-8, for which the max- 
imum value of A12 where state (0,1) can be stabilized 
is enhanced by ~ 45% compared to $m = case. We 
believe that these values can be considerably improved 
by exploring other configuration of the magnetic dot. 
However, the search of the optimum configuration is be- 
yond the scope of this work; we restrict ourselves here 
to the proof-of-concept. Although we present here only 



the (0, 1) case, we verified that the same conclusions ap- 
ply to other fractional states, but of course with different 
optimum magnetization values. 



V. CONCLUSIONS 

In summary, we have analyzed the stability of frac- 
tional vortex states in a mesoscopic superconductor with 
respect to changing its microscopic material parameters 
and to the magnetization state of a magnetic dot placed 
nearby. By sweeping a broad range of values of the par- 
tial density of states and of the Fermi velocities for a 
given coupling matrix, we have determined the stabil- 
ity regions of many different fractional vortex states. We 
could identify two main groups of such solutions: the first 
one comprises those states where the dominating band, 
i.e. that one with the highest Tcj, has a lower total wind- 
ing number; the second group comprises states where 
the dominating band has the highest winding number. 
As expected, states of the first group have, in general, 
lower energy and are easier to find because the dominat- 
ing band has a natural tendency to be more refractory 
to vortices. In addition, they can be found in thermo- 
dynamic equilibrium for the case of weak Josephson cou- 
pling, as previously reported in the literature J^ii^ How- 
ever, the states in each of these groups can exchange 
roles smoothly by sweeping the partial density of states. 
We have also demonstrated that the healing length of a 
band at zero temperature, defined as = hvj/VQW, is 
closely connected to the size of a vortex in that band and 
thereby to the ability of the surface barrier to superheat 
or supercool fractional states in applied magnetic field. 

We have further investigated the properties of frac- 
tional vortex states in inhomogeneous magnetic field. In 
particular, we exposed the superconducting disk to the 
stray field of a magnetic dot, placed coaxially above the 
sample. We demonstrated that such a configuration is 
capable of augmenting the stability region of fractional 
states, at least when the dot magnetization is parallel to 
the external homogeneous field. Moreover, the maximum 
value of interband coupling for which a given fractional 
state is stable can be increased considerably and reaches 
an optimal value at a well defined magnetization of the 
dot. This is an important result because fractional states 
are usually possible only in materials with very weak in- 
terband coupling. In other words, the here proposed idea 
enables observation of fractional vortex states in a con- 
siderably wider range of multiband materials. 

Finally, we point out that, to date, experiments de- 
termining vortex structures in two-band superconduc- 
tors have been restricted to bulk single crystals. How- 
ever, mesoscopic systems present a much richer variety 
of vortex states, including non-composite fractional vor- 
tices, which represent a clear, direct evidence of com- 
peting length scales in multiband materials. We believe 
our findings contribute to the construction of a roadmap 
guiding further search of suitable materials for imaging 
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these exotic vortex structures. In this respect, further 
investigation on other geometries of the mesoscopic su- 
perconductor and the magnetic dot aiming at the stabi- 
hzation of fractional vortices in wider parameter ranges 
would be very much welcome. 
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